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Abstract 

We introduce a new algorithm denoted Dsc^ to isolate the real roots of a univariate square-free polynomial 
/ with integer coefficients. The algorithm iteratively subdivides an initial interval which is known to contain 
all real roots of /. The main novelty of our approach is that wc combine Descartes' Rule of Signs and 
Newton iteration. More precisely, instead of using a fixed subdivision strategy such as bisection in each 
iteration, a Newton step based on the number of sign variations for an actual interval is considered, and, 
only if the Newton step fails, we fall back to bisection. Following this approach, our analysis shows that, 
for most iterations, we can achieve quadratic convergence towards the real roots. In terms of complexity, 
our method induces a recursion tree of almost optimal size 0{n ■ log(nT)), where n denotes the degree of 
the polynomial and T the bitsize of its coefficients. The latter bound constitutes an improvement by a factor 
of T upon all existing subdivision methods for the task of isolating the real roots. In addition, we provide 
a bit complexity analysis showing that Dsc^ needs only O(m^t) bit operations ' to isolate all real roots of 
/. This matches the best bound known for this fundamental problem. However, in comparison to the much 
more involved algorithms by Pan and Schonhage (for the task of isolating all complex roots) which achieve 
the same bit complexity, DSC'^ focuses on real root isolation, is very easy to access and easy to implement. 

Key words: Root isolation, Subdivision Methods, Exact Computation, Newton's Method, Descartes' Rule 
of Signs, Bit Complexity, Asymptotically Fast Methods 



1. Introduction 

Finding the roots of a univariate polynomial / is considered as one of the most important 
tasks in computational algebra. This is justified by the fact that many problems from mathemat- 
ics, engineering, computer science, and the natural sciences can be reduced to solving a system 
of polynomial equations which in turn, by means of elimination techniques such as resultants 
or Grobner Bases, reduces to solving a polynomial equation in one variable. Hence, it is not 
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surprising that numerous approaches are dedicated to this fundamental problem. We mainly dis- 
tinguish between (1) numerical and exact methods, and (2) methods to find all complex roots and 
methods which are especially tuned to search for real roots. The numerical literature lists many 
algorithms, such as Newton iteration or the Weierstrass-Durand-Kerner method, that are widely 
used and effective in practice but lack a guarantee on the global behavior (cf. [24] for discussion). 
In particular, the convergence and/or the complexity of the Weierstrass-Durand-Kemer method 
is still open. 

The work of A. Schonhage [29] from 1982 marks the beginning of the complexity-theoretic 
approaches. It combines a newly introduced concept denoted splitting circle method with tech- 
niques from numerical analysis (Newton iteration, Graeffe's method, discrete Fourier transforms) 
and fast algorithms for polynomial and integer multiplication. For the benchmark problem of iso- 
lating all complex roots of a polynomial / of degree n with integer coefficients of modulus 2^ 
or less, the proposed method achieves the record bound of d{r?T) bit operations. V.Pan and 
others [23,24] gave theoretical improvements in the sense of achieving record bounds simulta- 
neously in both bit complexity and arithmetic complexity, but the initial bound 0{n^x) on the 
number of bit operations has still remained intact. Common to all asymptotically fast algorithms 
is (as the authors themselves admit) that they are rather involved and very difficult to implement. 
The latter is also due to the fact that, in order to control the precision errors in the considered nu- 
merical subroutines, one has to carefully work out many details of their implementation. Hence, 
it is not surprising that, despite their theoretical richness, the asymptotically fast algorithms have 
so far not been used, or not proven to be efficient in practice; see [12] for an implementation of 
the splitting circle method within the Computer Algebra system Pari/GP. A further reason might 
be that the benchmark problem is inappropriate for most applications. For instance, in ray shoot- 
ing in computer graphics, we are only interested in the first positive root or in the real roots in 
some specified neighborhood. 

In parallel to the development of purely numerical methods, there is a steady ongoing re- 
search on exact subdivision algorithms, such as the Descartes method [7,9,15,18,26], the Bolzano 
method [4,5,6,17,33], Sturm Sequences [8,16,25] or the continued fraction method [2,30,31]. 
These methods from the exact computation literature are widely used in various algebraic ap- 
plications (e.g., cyUndrical algebraic decomposition), and many of them have been integrated 
into computer algebra systems (e.g.. Maple, Mathematica, Sage, etc.). In addition, their 
computational complexity has been well-studied [8,10,30,31,33], and many experiments have 
shown their practical evidence [13,14,26]. Current experimental data shows that a version of 
the Descartes method (i.e., the univariate solver in Rs based on [26], integrated into Maple) 
which uses approximate computation performs best for most polynomials, whereas, for harder 
instances, the continued fraction approach seems to be more efficient. With respect to the bench- 
mark problem, all of the above mentioned algorithms demand for d{n^X^) bit operations to 
isolate all real roots, hence they tend to lag behind the asymptotically fast algorithms by a factor 
of nx. Recently, it has been shown [28] that the bound on the bit complexity for the Descartes 
method can be lowered to d{ir'x^) when replacing exact computation by approximate compu- 
tation (without abstaining from correctness). This result partially explains the success of such 
a modified Descartes method in practice. However, as long as we restrict to the bisection strat- 
egy, it seems that the latter bound is optimal. We remark that Schonhage already made a similar 
observation: In the introduction of [29], he argued that "a factor inevitably occurs if nothing 
better than linear convergence is achieved". 
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In this paper, we introduce an exact and complete subdivision algorithm denoted Dsc^ to 
isolate all real roots of a square-free polynomial with integer coefficients. Similar to the classical 
Descartes method, we use Descartes' Rule of Signs to determine an upper bound v/ = var(/, /) 
for the number of real roots of the polynomial / within an interval / that is actually processed. 
However, instead of splitting / into two equally sized subintervals in each iteration, we consider 
a subdivision strategy that is based on Newton iteration: The analysis of the classical approach 
which exclusively uses bisection shows that the induced recursion tree is large (« nT) if and 
only if there exists a long sequences /i D /i D • • • D of intervals in the subdivision process, 
where v = v/j = • • • = v/^. Such a sequence implies the existence of a cluster of v nearby roots, 
and vice versa (cf Theorem 1). Hence, it seems reasonable to obtain a good approximation (i.e., 
an interval /' C / close to of such a cluster by considering a corresponding Newton step 
to approximate a v-fold root. Combining Descartes' Rule of Signs and a subdivision technique 
similar to the one as proposed by J.Abbott for the QIR method [1], we formulate a method to 
determine whether the so-obtained approximation /' should be kept or not. In case of success, 
we proceed with the considerably smaller interval /', whereas we fall back to bisection if the 
Newton step fails. Our analysis shows that, following this approach, we achieve quadratic con- 
vergence in most iterations. As a consequence, the induced recursion tree has almost optimal 
size 0(nlog(nT)) which improves upon the bisection strategy by a factor of T. We further pro- 
vide a detailed bit complexity analysis which yields the bound (5(n^T) for Dsc^. This matches 
the record bound achieved by the aforementioned asymptotically fast algorithms. 

We consider our contribution of great importance because of the following reasons: (1) Al- 
though the proposed method is rather simple, it achieves the best bounds known for the bit 
complexity of the problem of isolating the real roots of a polynomial. (2) In addition, it is much 
easier to access and also much easier to implement than the asymptotically fast algorithms that 
are available so far. In comparison to the existing practical methods for real root isolation, the 
modifications are moderate, and thus we expect that a careful implementation of our new ap- 
proach will outperform the existing ones. (3) Finally, our method can be applied to search for 
the real roots in some specified neighborhood of interest, a property which is not fulfilled by the 
algorithms as proposed by Pan and Schonhage. 

2. Overall Idea 

In this section, we first provide a high-level description of our new algorithm, and then outline 
the argument why this approach improves upon existing methods such as the classical Descartes 
method. For the exact definition of the algorithm and a detailed bit complexity analysis, we refer 
the interested reader to Section 3. Throughout the following considerations, let 

n 

fix) := £ Uix' e Z[x], with \ai\ < 2^ (2.1) 

(=0 

be a square-free polynomial of degree n with integer coefficients of bit-length t or less. We 
further denote zi, • • • ,Zn the complex roots of /, (j{zi) := min^-^,- Iz,- — z^l the separation ofzu and 
Of := min; C7(z;) the separation off. According to Cauchy's bound, the modulus of each root z, 
is bounded by 1 + 2^ < 2^+\ and thus, for the task of isolating the real roots of /, we can restrict 
our search to the interval J^o := (-2^+\2^+i). 

We consider an arbitrary root isolation method denoted ISO which recursively performs sub- 
division on /o in order to determine isolating intervals for the real roots of /. ISO uses a counting 
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Fig. 2.1. For any k with Q<k<n, tiie Obreshkojf discs Ci and Cf^ for / = {a,b) liave tlie endpoints of / on 
tlieir boundaries; their centers see the line segment {a,b) under the angle 2% / (k~\-2). The Obreshkojf lens 
Li; is the interior of Q nC;^, and the Obreshkojf area A/; is the interior of Q UQ^. Any point (except for a 
and b) on the boundary of Aj. sees / under an angle n/(k + 2), and any point (except for a and b) on the 
boundary of Lj. sees / under the angle 7t — 7t/{k + 2). We have L„ C . . . C Lj C Lq and Aq <Z A\ <Z . . . <Z A„. 
The cases k — and = 1 are of special interest: The circles Co and Cq coincide. They have their centers 
at the midpoint of /. The circles Ci and Cj are the circumcircles of the two equilateral triangles having / as 
one of their edges. We call Aq and Aj the one and two-circle regions for /. 

function var(/,/) for m, the ntunber of roots within an interval / C J^o. where v := var(/,/) e N 
fulfills the following properties: 

(PI) V is an upper bound for m (i.e., v > m), and 
(P2) V has the same parity as m (i.e., v = m mod 2). 

The latter two properties imply that var(/, /) = m if v < 1 . Hence, in each step of the recursion, an 
interval /is stored as isolating if v = 1, and / is discarded if v = 0. If v > 1, / is subdivided 

(according to some subdivision strategy) into subintervals /i = {a,Xi), h = (/li,/l2), . . . ,// := 
(A/_i ,£)), where 1 < / < /q, /o € N is a global constant, and the A, are rational values. In order to 
detect roots at the subdivision points A,, we also check whether /(A;) = or not. Tiso denotes the 
recursion tree induced by ISO. Throughout the following considerations, we often treat nodes of 
ISO and intervals produced by ISO as interchangeable. 

The definition of the counting function var(/,/) is based on Descartes' Rule of Signs: For an 
arbitrary polynomial p = Y!i=o Pi^' & ^[x], the number m of positive real roots of p is bounded 
by the number v of sign variations in its coefficient sequence (/?o, • . . and, in addition, 
v = m mod 2. In order to extend the latter rule to arbitrary intervals I = {a,b), the Mobius trans- 
formation X I-)- which maps (0, +«>) one-to-one onto / is considered. Thus, for 



and var(/,/) defined as the number of sign variations in the coefficient sequence (cq, . . . , c„) of //, 
var(/,/) fulfills the properties (PI) and (P2). Because of the latter two properties and the fact that 
we never discard intervals that contain a real root of /, correctness of ISO follows immediately. 

The classical Descartes method (Dsc for short) is a subdivision method which uses bisec- 
tion in each iteration, that is, in each step, we have I =Iq = 2, h = (a,Ai) := {a,m{I)) and 




(2.2) 
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h = (■^1 , b) = {m{I) , b), with m{I) := {a + b)/2 the midpoint of /. Termination and complexity 
analysis of Dsc rest on the following theorem: 

Theorem 1 ([19,21]). Let I = {a,b) be an open interval and v — var(/,/). If the Obreshkojf 
lens Ln-k C C (see Figure 2.1 for the definition ofL^-k) contains at least k roots (counted with 
multiplicity) of f, then v>k. If the Obreshkojf area A^dC contains at most k roots (counted 
with multiplicity) ofk, then v<k. In particular, 

fP3) # of roots off in L„ < var(p,7) < # of roots off in A„. 

We remark that the special cases k = Q and k= \ appear as the one- and two-circle theorems 
in the literature [3,9,15,20,22]. For the Descartes method. Theorem 1 impUes that no interval / 
of length w{I) < (7/ is split. Namely, its one-circle region Aq cannot contain two or more roots. 
If Aq contains no root, then vai{f,I) = 0. Otherwise, Aq contains one real root, and thus the two- 
circle region Ai contains no non-real root. Hence, in the latter situation, we have var(/,/) = 1 
by Theorem 1 . We conclude that the depth of the recursion tree Tbsc induced by the Descartes 
method is bounded by logw(j^o) +logaJ^ = T-|-logoy ^ +2. Furthermore, it holds (see [9, 
Corollary 2.27] for a self-contained proof): 

Theorem 2. Let I be an interval and I\ and h be two disjoint subintervals of I. Then, 
(P4) var(/,/i) +var{f,h) < yw{f,I). 

According to Theorem 2, there cannot be more than n/2 intervals / with var(/,/) > 2 at 
any level of the recursion. Hence, the size of Tbsc is bounded by n(T + logO'^' +2). Using 
Davenport-Mahler bound, one can further show [9,28] that logcT^' = 0{n{\ogn + t)), and thus 

the bound for | Tbsc | writes as d{n^x) . A more refined argument [9] yields | Tdsc | = Oinx) which 
is optimal for the bisection strategy. 

In the next step, we study the situation where the recursion tree Tdsc for the classical Descartes 
method is large. We then introduce our new algorithm which we denote Dsc^ due to the quadratic 
convergence in most steps. DSC^ is a variant of the Descartes method which adaptively addresses 
the latter situation via combining Newton iteration and bisection. We also sketch the argument 
why this approach improves upon Dsc. The following definition is essential for the argument; 
see also Figure 2.2: 

Definition 1. Let T be a node (interval) in the recursion tree Tiso induced by some subdivision 
algorithm ISO. We call T terminal if var(/,T) < 1 . A non-terminal interval I with children I\,...,Ii 
is called special HI = J'q (i.e., I is the root of Tiso), or 

^^AfJj) < var(/,T) for all j = 1, . . . 

According to (P4) in Theorem 2, we have var(/,Tj) < var(/,T) for each I. Thus, a 
non-terminal node I different from J^o is non-special if and only if, for one of its children, we 
count the same number of sign variations as for I and, for all other children, we count no sign 
variation. In total, there exist n' special nodes, where n' < var(/, J^o) < n. Namely, when we 
subdivide a special interval which is not the root of the recursion tree, the non-negative value 
ji := L/ var(/,T) — #{T : var(/,T) > 0} decreases by at least one, where we sum over all leafs in 
the actual iteration, and is initially set to ;U = var(/, J^o) — 1 <n—l. We denote the special 



5 



Fig. 2.2. The left figure shows the subdivision tree Ibsc induced by the Descartes method, where, for each 
node /, the number var(/,/) of sign variations is given (e.g., var(/,yi) = 9, or var(/,y4) = 2). The colors 
red, black and blue indicate the (1) special nodes, (2) terminal nodes, and (3) all other nodes (non-special 
and non-terminal), respectively. The right figure shows the subtree T^^^ obtained by removing all terminal 
nodes. The non-special nodes in T^^q partition into maximal connected components T4 and Te,. 

nodes by 7i, . . . and assume, w.l.o.g., that vv(7,) > w{Jt) if i < k. In particular, J\ = J^o- We 
define T^^^ the subtree of Tiso obtained from Tiso via removing all terminal nodes. Then, Tj*,-, 
partitions into 

(1) the special nodes Ji,...,J„i (red dots in Figure 2.2), and 

(2) subtrees T, C T^^q, with i — 2, . . . ,n' , consisting of all non-special nodes / G T^^^ with 
Jj C / and <f. I for all special nodes 7^ with 7^ D 7,- (blue dots). 

From our definition of a special node, it follows that each 7] constitutes a chain of intervals 
h Z) ■ ■ ■ D Is that connects two special nodes. More precisely, 7; connects J, with J^, where is 
the special node of minimal width that contains 7, . Since each interval has at most Iq children, 
I Tiso I is bounded by {Iq + I) ■ \ T^^q\ - Hence, we have 

o(|r,so|) = o(|r,;ol) = o(n' + |:|7;-|) = o(«) + o(|:|7;.|). (2.3) 

(=2 i=2 

The latter consideration shows that the size of the subdivision tree mainly depends on the 
length of the chains T,. For the Descartes method, it might happen that some of these chains are 
very large (i.e., | « nx) which is due to the following situation (see also Figure 2.3): For a 
polynomial / as in (2.1), it is possible that there exists a ^ e M and a very small, complex neigh- 
borhood (of size e « 2^"^) of ^ that contains a cluster ^ of v nearby roots of /. Thus, separating 
these roots from each other via bisection requires at least loge^' « nT steps. Furthermore, due 
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Fig. 2.3. The long chain /i D /2 D • • • D /j of intervals in T{5) with var(/,/i ) = •••= vai{f,Is) = 2 corre- 
sponds to a large number of bisection steps to isolate two very nearby roots from each other; see the figure 
on the right with the graph of / over /] . 

to (P3) in Theorem 1, there exists a long sequence /i D /2 D • • • D of non-special intervals with 
^ G Ij for all j, and thus the number 

V var(/,/i) = var(/,/2) = • • • = var(/,/,) 

of sign variations does not change for the intervals in this sequence. Namely, for each Ij in the 
above sequence, the Obreshkoff lens L„ contains Vice versa, according to Theorem 1, such a 
long sequence of non-special intervals implies the existence of a cluster consisting of v nearby 
roots as above because the Obreshkoff area A„ of each Ij must contain at least v roots Since a 
cluster of v nearby roots at ^ behaves very similar to a v-fold root at ^ , it seems reasonable to 
obtain a good approximation of by considering Newton iteration instead of bisection. Namely, 
for a polynomial p{x) E M.[x] with a v-fold root at ^ and a starting value xq sufficiently close to 
it is well-known from numerical analysis that the sequence recursively defined by 

converges quadratically to ^ . Unfortunately, when isolating the roots of /, the situation differs 
considerably from the latter one: First, the above result only holds for a v-fold root ^ and does 
not directly extend to a cluster of v roots near £, . Second, in an early stage of the subdivision 
process, the existence of such a cluster 'yf is not guaranteed, and even if one exists, we do not 
know what "sufficiently close to t^" means in this situation. 



^ The thoughtful reader may notice that the latter two statements are not completely rigorous: In particular, in the special 
case where ^, and thus also the cluster is very close to one of the endpoints of some Ij, it might happen that some of 
the roots are not considered by the counting function var(/,/j) since they are located outside the Obreshkoff lens/area. 
We will address this issue in our algorithm as defined in Section 3. 
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In order to address the above mentioned problems and to finally turn the purely numerical 
Newton method into an exact and complete algorithm, we propose the following approach: Let 
V = var(/, /) be the number of sign variation for an actual interval / = (a, ^) in a certain iteration. 
Then, we consider this as an indicator that there might exist a cluster of v nearby roots. Thus, 
we compute A := t ~ v ■ f{t) / f [t] for some t G [a,b] (e.g., an endpoint of /) and consider an 
interval I' = {a' ,b') Cl of width w(/') <C w(/) that contains A. If var (/,/') = v as well, we keep 
/' and discard the intervals {a, a'] and [b',b). Otherwise, we split / into two equally sized intervals 
Ii := {a,m(I)) and h '■= {m{I),b) and finally check whether f(m(I)) = or not. Following this 
approach, no root is lost and intervals are at least bisected in each iteration. Furthermore, if a 
cluster of nearby roots actually exists, we can hope to achieve fast convergence to this cluster 
when choosing /' in an appropriate manner. In our algorithm, we choose /' in a similar way as 
proposed by Abbott [1] for the task of further refining intervals which are already isolating for 
an ordinary root. Namely, we decompose / into a certain number Nj of subintervals and pick the 
subinterval /' of size w{I) /Ni which contains A. If var(/,7') = v, then we keep I' and decompose 
/' into N[i = Nf subintervals in the next iteration. Otherwise, we continue with the intervals 
/i = {a,m{I)) and/2 = {m{I),b) which are now decomposed into only A'/j =Nj^ :=max(4, 
many subintervals, etc. 

In the next section, we give the exact definition of our new algorithm, and we show that it 
induces a subdivision tree of considerably smaller size than Tdsc- In particular, it turns out that 
the size of each 7^ c T^^^i is bounded by 0(log« + logT) which is due to the fact that, for most 
iterations, we have quadratic convergence to the real roots, and the width of each interval is lower 
bounded by 2^''("^) ; see Lemma 1 and Theorem 4 for proofs. Hence, according to (2.3), the size 
of the overall recursion tree is bounded by 

0{n' ■ (lognr)) = 0(var(/, .^0) ■ (lognr)) = 0{n ■ (lognr)). (2.4) 

The latter result particularly shows that the size of the recursion tree is directly correlated to the 
number n* of non-zero coefficients of / because instead of considering J^o = (— 2^+^ , 2^+' ), we 
can start with (—2^+^ , 0) and (0, 2^+^ )), and the total number of sign variations counted for both 
intervals is upper bounded by 2 • n*. 

For the bit complexity of our algorithm, we have to consider the costs for computing the 
polynomials fj{x) = (1 +x)" ■ f{{ax + b)/{\ +x)) as defined in (2.2), where / = {a,b) is an 
interval to be processed. In Section 3.3, we will show that the costs for the latter step are domi- 
nated (up to constant factors) by the computation of /(x + a). For / G Ti, the endpoints of / are 
dyadic numbers of bitsize 0(T + logw(7,)"^) or less, and thus the computation of // demands 
for d{rt}{T + \ogw{Ji)~^)) bit operations. In Lemma 6, we prove that we can order the roots 
zi, . . . ,z„ in a way such that logw(7,)"' < logfj(z,)~^ -l-O(log'^n) for all j = 1, . . . ,n'. Then, it 
follows that computing // demands for at most 

d{n^{T + \ogw{Ji)-^)) = d{n^{T + \oga{zi)-^+\og^n)) = d{r?{x + \ogC!{zi)-^)) 

bit operations. Thus, for the total cost, we obtain the bound 

d{Tv'T + n^Y,^oga{ziT^) = O(n^T) 

1=1 

since Y4=\ logc7(z/)"^ = 0(nx{\-Ognx)) = d{nx) according to Lemma 19 in [28]. 
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3. Algorithm and Analysis 

3.1. The Algorithm 

We first present our new algorithm denoted Dsc^. For pseudo-code, we refer to the Appendix. 

DSC^ maintains a list of active intervals / with corresponding integers Nj = TF' , ni € N> i , 
and a Ust of isolating intervals, where we initially set := {{J'o^A)} := {((-2^+\2^+i),4)} 
and ^ := 0. For (/, A^/) € / = (a, b), we proceed as follows: We remove / from ^ and compute 
the number v :— var(/,/) of sign variations for / on /. 

(1) If V = 0, we do nothing (i.e., / is discarded). 

(2) If V = 1, then I isolates a real root of /. Thus, we add / to the list of isolating intervals. 

(3) For V > 1, we proceed as follows: 

(a) For J = 1,2, let 

Bi:=(a,a+'^) and 82:= {b-'^^b) (3.1) 

be the left- and rightmost interval of size w{I)/Ni contained in /, respectively. We 
compute V, := vaT{f,Bi): If one of the values vi or V2 equals v, then the corresponding 
interval 5, (at most one of the two disjoint intervals fi, fulfills var(/,fi() = v) contains 
all roots of / within /. Hence, we keep /' := B,, discard /\/', and set Np :=Nf. That 
is, (I'.N,:) := (Bi.Nj) is added to s^. 

(b) If both values vi and V2 differ from v, we compute 

/(«) f{b) 
If I contains a cluster ^ of v nearby roots and a (or b) has "reasonable" distance to 
'i', then Xi (or A2) constitutes a considerably better approximation of than a (or b). 
We check whether this is actually the case: For ; = 1,2, we first compute the point 
a+kf ^ , with fe,- e {2, . . . , 4Nr — 2}, which is closest to A,- (if there exist two equally 
close points, we choose the one with smaller index). In more mathematical terms, 

ki :=min(max([4A^/(A, -fl)J,2),4A^/-2). (3.3) 

Then, we define /• to be the interval of length w{I)/Ni centered at the subdivision 

point a + ki- w{I) /4Ni, that is, 

4 : = (a -f - 2) . ^ , a + (^,- + 2) . ^ ) C /. (3 .4) 

In particular, if a + w{I)/4Nj <Xi<b — w{I) /4-Nj, then /' contains A,-, and A,- has dis- 
tance at least \v{I)/4Ni to both endpoints of /■. Now, we compute v\ := var(/,/J ) and 
Vj := var(/,/2). If one of the two values v\ or v'2 equals v, we keep the corresponding 
interval /' := /• with v' = v (if we count v sign variations for both intervals I\ and I2, 
we just keep I[) and discard /\/'. Finally, we add {I',Njt) := {I',Nf) to £/. 

(c) If all values v\, V2, v[ and differ from v, then we consider this an indicator that 
there is either no cluster of v nearby roots or that such a cluster is not separated well 
enough from the remaining roots. Hence, in this situation, we fall back to bisection. 
That is, we split / into two equally sized intervals /i = (a,m{I)) and I2 = {m{I),b) 
and add (/i,max(4, VA^)) and (/2,max(4, v^)) to ^. Finally, if /(m(/)) = 0, we 
also add [m{I),m{I)] to 
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Correctness and termination of DSC^ are obvious because our starting interval J^o contains 
all real roots of /, we never discard intervals (or endpoints) that contain a root of /, and intervals 
are at least bisected in each iteration. In addition, we obtain the following bounds on the width 
of the intervals / and the corresponding numbers Nj produced by Dsc^: 

Lemma 1. For each interval / produced by Dsc^, we have 

2^+2 > w{I) > aj ■ 2-2(^+2) = 2-^(«^) and 4 < AT/ < 2^(^+2) • af = 2^^"^l 

In particular, for I GTj (see Section 2 for the definition of the subtree Tj c Tj^gQi), we have 

2^+2 > w(/) > w{Ji) and 4 < AT/ < 22(^+2) • w(y^)"^ 

with /, the special node corresponding to 7] . 

Proof. The inequalities 2^+2 > ^(^/^ ^nd Nj > 4 arc trivial. For Nj > 4, there must exist an 
interval J Dl with Nj = \/Ni, and J was replaced by an interval 7' 3 7 of size w{J) /Nj. Since J 
is non-terminal, J' is also non-terminal because var(/,7') = var(/,7) > 1. Thus, oy < w{J') = 
< 2^+2/v^. This shows the upper bound for Ni. For the lower bound for w(/), we 
consider the parent interval J of I. Since J is non-terminal, we have w{J) > oy, and thus w{I) > 
w{J) /Nj >af aj- 2-2(^+2) . For / G 7]-, the bounds for w(7) are trivial, and, in completely similar 
manner as above, we conclude that 2''+^/^jNi > w{Ji). □ 

Throughout the following considerations, we call a subdivision step from 7 to 7' C / quadratic 
if ^(7') — w{I)/N/, and we call a subdivision step linear if 7 is split into two equally sized 
intervals 7i and 72. In a quadratic step, the integer Ni is squared whereas, in a linear step, Np := 
max (4, for each of the subintervals 7' =Ii/2- 

3.2. Analysis of the Recursion Tree 

In this section, we prove that the size of each of the subtrees 7} C T^sc^ defined in Section 2 
is bounded by 0(logn + log t). We first have to investigate into the following technical lemmata: 

Lemma 2. Let w, vJ € K+ be two positive reals with w > W , and let m G N>i be a positive 
integer. The sequence (s,),gN>i '•= ((•^i)i())(6N>i is recursively defined as follows: s\ = (x\,n\) = 
(w,ot), and 

(^fi,max(l,n,-_i-l)), if|^<w', 

where A^, := 22"' and i>2. Then, the smallest index j'o with x,^ < w' is upper boimded by 8(ni -|- 
loglogmax(4,^)). 

Proof. Throughout the following consideration, we call an index / strong (S) if x,/A^, > w' and 
weak (W), otherwise. If w/4 < w', then each ; > 1 is weak, and thus i'o < 3. For w/4 > w', let 
k e N>i be the unique integer with 

2-2^^'<w7w<2-2\ 

Then, k < log log ^, and since x, < x,-i/2 for all i, there must exist an index i which is weak. 
Let kf denote the smallest weak index. 



{xi,ni) := 
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Claim 1: k' <k+l 

Assume otherwise, then the indices 1 to are all strong. Hence, 

and n^+i > 1. It follows that A: + 1 is weak, a contradiction. 

Let us now consider the subsequence = kf,kf + l,...,io — 3: 
Claim 2: y contains no subsequence of type ...SS... or ...SWSWS... 

If there exists a weak index / and two strong indices /+ 1 and i + 2, then Xi/Ni > Xi+2/Ni > 
Xi+2/Ni+2 > contradicting the fact thatx,/iV,- < w'. Since y starts with a weak index, the first 
part of our claim follows. For the second part, assume that i, i + 2 and ? '+ 4 are strong. Then, i + 1 
and J ' + 3 are weak, and thus 

w' < Xi+4/Ni+4 < Xi+2/{Ni+2-Ni+A) < Xi/{Ni-Ni+2-Ni+4) = Xi/Nf = Xi+i/Ni+i 

contradicting the fact that j ' + 1 is weak. 
Claim 3: If i is weak and i < io, then n,- > 2. 

Namely, if / is weak and n,- = 1, then Xi/4 = xi/Ni < W , and thus x,g_i < w' which contradicts 
the definition of jq. 

We now partition the sequence ,5^ into maximal subsequences , ■ ■ - T^r such that each 
S^j, j = l,...,r, contains no two consecutive weak elements. Then, according to our above re- 
sults, each S^j, with j < r, is of type W, WSW, or WSWSW. The last subsequence S^r (with last 
index /q - 3) is of type W, WS, WSW, WSWS, or WSWSW. After each Sj, with j < r, the num- 
ber rii decreases by one, and thus we must have r <ni+k' since we start with rn^i = m + k' — I 
and, in addition, n, > 2 for all weak i. Since the length of each is bounded by 5, it follows 
that 

iQ = iQ-3 + 3 <k' +5r + 3 <5{ni+k')+k' + 3 < 8{ni+k). 

□ 

Lemma 3. Let / = {a,b) an arbitrary interval, A„ the corresponding Obreshkoff area and L„ the 
Obreshkoff lens for /. Then: 

(1) For I' = {a' ,b') C I with a a' and b b', the Obreshkoff area for /' is completely 
contained within the lens L„ if 

miR{\a-a'\,\b-b'\) >Sn^w{l'). 
In the latter situation, it holds that 

\x-^\ > j^ - {nan{\a-a'\,\b-b'\)-Snw{l')) 

for all x^L„ and all ^ eA'^. 

(2) For/' = {a',b') with /'("!/ = 0, the Obreshkoff area for/' does not intersect A„ if 

dist(/,/') > 4«2 .min(w(/),w(/')), 
where dist(/,/') denotes the distance between the two intervals / and /'. 
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Fig. 3.1. On the left figure, c denotes the topmost point of the Obreshkoff lens L„ for / = {a^b). If 
\a — a'\ < w{I)/2, then the distance from a' to the boundary of L„ is bounded by the distance 5 from a' 

to etc. The radius / of the Obreshkoff discs Cj, and cj, for /' — {a' ,h') is bounded by n ■ vv(/') due to the 
extended Sine Theorem. The right figure shows the Obreshkoff areas for the intervals / and /', respectively. 

Proof. (1) In a first step, we compute the radius r' of the Obreshkoff discs Cj, and c'„ for the 
interval/' = {a',b'): A point ^ on the boundary of the Obreshkoff area A J, (except a' and b') sees 
/' under an angle 7= ;r/(n + 2); see Figure 2.1 and 3.1. Hence, from the extended Sine Theorem, 
it follows that 

win_ win ^ (n + 2)w(/) 

2sm(7) 2sm{7c/{n + 2) n ^ ' 

since sinx > x/2 for all x e (0,;r/4]. In particular, each point z within the Obreshkoff area Aj, 

has distance at most 2r' < 2« • w(/') from any point within /. W.l.o.g, we assume that |fl — a'| < 
\b — b'\. Then, the distance from a' to the boundary of the Obreshkoff lens L„ for / is bounded 
by the distance 5 from d to the line ac, where c denotes the topmost point of L„. Since ac 
intersects the x-axes in an angle of ;r/(2(n + 2)), we have 5 = ja — fl'| sin;r/(2(n + 2)) > \a — 
d\l{\n). Thus, a;, CL„ if |a-a'|/(4n) > 2n-w(/') or \a-cJ\ > 8n^w{I'). For \a-a'\ > \b-b'\, a 
similar argument shows that Aj, C L„ if \b — b'\ > Sn^w{I'). In addition, if min(|fl — a'|, \b — b'\) > 
8n^vv(/'), then each point ^ <eA'^ has distance at least 5 — 2nw(/') >rain{\a — a'\,\b — b'\)/{4n) — 
2nw{I') from any point located outside of L„. 

(2) W.l.o.g., we can assume that w(/') < w{I) and a' > b. Let L be the line passing through b 
which intersects the x-axes in an angle of ;r/(n + 2). Then, the upper part of the Obreshkoff area 
A„ lies completely on one side of this line. Now, if A^ lies completely on the other side of L, 
then, by symmetry, A„ and A'^ do not share a common point. We have already argued that Aj, is 
contained within the disc of radius 2nw{I') centered at a' . Hence, if the distance 5' := dist(fl',L) 
from a! to L is larger than 2nw(/'), then A„ n Aj, = 0. We have 5' = |a' - b\ sm{n/{n + 2)) > 
\a' — b\/2n = dist(/,/') /2n, and thus our claim follows. □ 

We now turn to the analysis of the subtrees Tj C Tdjc^ defined in Section 2. There, we 
have already argued that each 7^ constitutes a chain of intervals /i = D = (^25^2) 3 

• • • D Is = {as,bs) "connecting" the special node 7; with the special node Jk of minimal length 
that contains 7,. In the following Theorem, we will show that, for aU but (9(logn) many j, the 
sequence {w{Ij),nij) = (w(/y),loglogA^/^.) behaves similar to the sequence {xj,nj) as defined in 
Lemma 2. As a result, we obtain the following bound on |7]|: 

Theorem 4. For each special node 7„ the corresponding subtree 7} C Tj^gc2 has size 

|71| =0(logn + logT). 

Proof. We first consider the special case where ai = 02 = ■ ■ ■ = Os, that is, in each subdivision 
step, the leftmost interval has been chosen. Since v = var(/| ) = •••= var(/5). Theorem 2 imphes 
that var(/, /) = v for each interval / with C / C /i . In particular, if w{Ij) /Njj > w{Is), we count v 
sign variations for the interval fii = {aj,aj + w{Ij)/Ni.) = {as,as + w{Ij)/Ni.) as defined in (3.1). 
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Thus, the subdivision step from /, to Ij+i is quadratic in this case. Then, for j — 1, ... ,i — 1, the 
sequence {w{Ij),njj) coincides with the sequence {xj,nj) as defined in Lemma 2, where w := 
w(/i),w' := w{Is) andni = m := n/j . Namely, \fw{Ij)/Ni. > w' , wehaveH'(/y-|-i) =w{Ij)/Ni. and 
n/j.^j = 1 +ni., and, otherwise, we have w{Ij+\) = w{Ij)/2 and n/^^j = max(l,n/^. — 1). Hence, 
according to Lemma 2, it follows that s is bounded by 

8(n/j +loglogmax(4,w(/i)/w(/^))) = 0(logn + logT), 

where we used the bounds on «/, and w{Is) from Lemma 1. An analogous argument shows the 
same bound for s in the case where bi=b2 = -bs- 

We now turn to the more general case, where a\ ^ a, and b\ ^ bg'. Let e {1, . . . be the 
smallest index with a^j ai and bg^ ^ b\ . Then, due to the above argument, is bounded by 
0(logn + logT). Furthermore, min(|ai —as^\^\b\—bs^\) > >v(/jj)/4, and thus 

min(|ai -aj\,\bi-bj\) > 2^-'^-^!]) 

1 /2^-*i-4 \ 
^ —■{mm{\ai-aj\,\bi-bj\)-Snw{Ij))>w{Ij)l 8nj for all j>si. 

Hence, with^i := + [log(16n^)] +4 = +0(logn), this yields 

1 2 
— ■ {nun{\ai—aj\,\bi—bj\) — Snw{Ij))>Sn w(/,)forall j>S2- 

Then, from Theorem 1 and Lemma 3, we conclude that, for j >S2, the Obreshkoff area A^-'^ for 
Ij contains exactly v roots zi , . . . ,Zv of / because the Obreshkoff lens L^^' for /j contains at most 
V roots and aI^^ C lI^'' contains at least v roots. In particular, the Obreshkoff area aI,^^ for must 

(A 

contain zi, . . . ,Zv. In the proof of Lemma 3, we have already argued that each point within A„ ' 
has distance less than 2nw{Is) from any point within Ig, and thus 

|x — z,| < 2nw{Is), for all i=l,.. . ,v and allx e Ig- (3.5) 

The remaining roots Zv+i,-- ■ ,z« of / are located outside the Obreshkoff lens Ln ^ for /i, and 
thus their distance to an arbitrary point within Ij is larger than 8n^w(/s2). Namely, according 
to Lemma 3, the distance from any of the roots Zv+i, ■■■■Z„ to an arbitrary point within I^^ is 
lower bounded by ^n^w{Is2) and contains Ij. The following consideration further shows the 
existence of an 53 = S2 + 0{logn + logT) such that w(/y) < w{Is2)/Nij for all 7 > S3, and thus 

\x-Zi\ >8n^A^/^w(/y)foralli = v+l,...,n, j>S3, andallxS/;-: (3.6) 

Due to Lemma 1, we have A'/^.<A?max:= r22(^+2)/oy] =2<5Mforall 7. Thus, if the sequence/s^, 
Is2+\ starts with more than mmax := loglogA^^ax + 1 = 0(logM+logT) consecutive linear sub- 
division steps, thenNj., =4andw(/y/) <w{Is2)/4 = w{Isi)/Ni., for some /< S2 + log log mmax- 
Otherwise, there exists a / with S2 < / < S2 +wimax such that the step from Iji to Iji^i is 
quadratic. Since the length of a sequence of consecutive quadratic subdivision steps is also 
bounded by mmax. there must exist a /' with / + 1 < /' < / -|-OTmax + 1 such that the step from 
//'_! to If is quadratic, whereas the step from /y// to/y//_|_i is linear. Then, A'/^.,,^^ = ^Ni.„ =Ni.ii_^ 
and 

w(/,»+i) = w{If)/2 = w{Ij„.,)/{2Ni.„J < w{Ig2)/Ni.„^,. 
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Hence, in both cases, we have shown that there exists an sj, < S2 + 2wmax + 1 = 0(logn + log t) 
with w{Isj) < w{Is^)/N[^^. Then, by induction, it follows that w{Ij) < w{Is^)/Nj. for all j > 
which shows (3.6). 

We are now ready to show that the subdivision step from Ij to is quadratic if j > sj, 
and w{Ij) > 6SnNijw{Is): Namely, if the latter two inequahties hold, then one of the endpoints 
(w.l.o.g., we assume this point to be aj) of Ij has distance at least w{Ij)/2 > 34nNijW{Is) from 
as- Thus, the distance from aj to any of the roots zi , . . . ,Zv is larger than 34nNijW{Is) — 2nw{Is) > 
32nNijw{Is). In addition, we have \aj — Zi\ > Sn^Nijw{Ij) for all « > v + 1 due to (3.6). Thus, 
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where we used that f'{a)/ f{a) = Y!i=i{'^ — Zi)~^ for all a e C with /(a) ^ 0. In completely 
analogous manner, we show that 



{aj-as)f'{aj) 



> 1 



1 

Wj- 



This yields the existence of an e e : 



. with |e| < l/{SNj) < 1/32 and i • ("^ ^"^j-^jK' = 1 +e. We 
can now derive the following bound on the distance between the approximation obtained by the 
Newton iteration and a.: 
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(3.7) 



If as > bj - w{Ij)/Nij, then (a,,/?,) C B2 = {bj - w{Ij)/Ni^,bj) (cf. Step 3 (a), (3.1) in our 
algorithm for the definition of B2), and thus var(/,B2) = v. Hence, in this case, we keep Ij+i = B2 
which has size w{Ij)/Nij. If < bj — w{Ij)/Nij, then according to (3.7) we have 



V ■ . e [as 



INi, 



AN,. 



)• 



It follows that the interval I[ as defined in Step 3 (b), (3.4) of our algorithm contains (as,bs), and 
thus var(/,/J) = v. This shows that the subdivision step from Ij to Z^+i is quadratic. 

We now consider the sequence {w{Is^+i),ni^^^^), for i = where i' is defined as the 

largest index with w{Is^+ii) > 68nw(/v). Then, our above argument implies that the sequence 
(w(/i3+/),n/j^_^,.)i<,<,'^i coincides with the sequence (x,-,n,)i<,<,/_i as defined in Lemma 2, 
where m = m = ni^^_^^ and w' := 68nw(/j). Namely, if w{Isj+i)/Ni^^^. > w', then 68nw(/s) < 

!+«/„+;. whereas, for w{Is^+i)/Ni . < w', we have 



w(/i3+,+i) = w{Is3+i)/Nsr^+i and n/^ 



loglogmax(4,w(/s3+i)/w')) = 0(logn+ log t). Hence, thereexistsans4 = s3+j' + l = 0(log«+ 



max(«/^ . — 1, 1). It follows that is bounded by 8(ni + 



14 



logr) with w{Ij) < 68nw(/v) for all j > S4. Finally, this shows that s is upper bounded by 

S4 + log(68n) = 0(logn + logT). □ 

Combining the latter theorem and (2.3) immediately yields the following result on the size of 
the induced recursion tree: 

Theorem 5. For a polynomial / of degree n with integer coefficients of bitsize T, the algorithm 
Dsc^ induces a recursion tree T-^^^i of size 

I^Dsc2 1 = var(/,/o) • 0(logn + log t) = 0{n ■ (logw + logr)), 
where Iq (— 2^+^ , 2^+^ ) denotes the initial interval known to contain all real roots of /. 

3.3. Bit Complexity Analysis 

We now derive an upper bound for the number of bit operations that are needed to determine 
isolating intervals for the real roots of /. We will show that, in each iteration, the costs are dom- 
inated (up to a constant factor) by the computation of the polynomial fi — (x+l)" ■ f{{ax + 
b)/(x+l)) as defined in (2.2). The costs for this step mainly depend on the bitsize of the end- 
points, and thus on the length of the interval / = {a, b). The following Lemma provides a lower 
bound on the width of the special nodes 7,-, and thus also for the nodes / e 7], in terms of the 
separations of the roots zi , . . . , z„: 

Lemma 6. For a polynomial / as defined in (2.1), we can order the roots of zi, . . . ,z„ of / in 
way such that 

^(y.)>^.„-5-2iog«^ 

where « = !,...,«' and Ji,...,J„' denote the special nodes in the recursion tree T^^^i . In addition, 
the endpoints of an arbitrary interval I £T (i) U {7, } are dyadic numbers which can be represented 
by 0(T-|-log(7ri -|-log^n) bits. 

Proof. We first order the roots of / with respect to their separations, that is, > ••• > (T(zn). 

For a fixed G E IR+, let k be defined such that (y{z„-k+i) < <^ ^ <^{z„-k), that is, exactly the 
k roots Zn-k+iT ■ ■ fZn have separation less than o. We further denote /i ;= 7,j ,...,/« := Ji„ the 
special nodes such that 

) < w^in : = ^ • n-5-2i°g« for alH = 1 , . . . , m, 

and each special node 7, G ^dsc^ which contains // has width w(y,) > Wmin- In addition, we denote 
V/ := var(/,//) > 2 the number of sign variations that we count for Since the intervals Ii are 

disjoint, we have vi H h < n, and thus m < n/2. According to Theorem 1, v; is a lower 

bound for the number of roots within the Obreshkoff area A^'' for Furthermore, in Lemma 3, 
we have proven that any two points within a['' have distance less than 4nw(//) < a, and thus 
each root contained in An^ must be one of the k roots Zn-k+ i,---,Zn. Let Si be the set of all roots 
which are contained in A^J^ . 

Let us first consider the case, where the Obreshkoff areas A^^ are pairwise disjoint. Then, the 
subsets Si C An ^ are also pairwise disjoint, and thus 
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We now turn to the more general situation, where some of the A„ ' may overlap. If this happens, 
then it is possible that, in total, some of the roots contained in these areas are counted more than 
once. Hence, we propose the following approach: In a first step (see the subsequent construc- 
tion), starting with the intervals /i , . . . ,/,„, we iteratively merge intervals whose corresponding 
Obreshkoff areas overlap until we finally obtain intervals such that 

• the intervals /i ,...,/„ are covered by /j , . . . , , 

• each /; has width w(//) < w^in • and 

• the Obreshkoff areas A„ for 7^' are pairwise disjoint. 

Then, using the same argument as above yields vi H h Vm < v' := Y!iLi var(/,/^') < k, where 

the latter inequality follows from the fact that the Obreshkoff area A„ for each 7^' contains only 
roots with separation less than 4nw{I'i) < 4nWmin • n'^+^'^s" < a{zn-k)- It remains to show how 
to construct the intervals /;': We start with a list of active intervals, where we initially set 
^ := In each iteration, we pick two intervals I = {a,b) and 7' = {a\b') from £/ 

whose corresponding Obreshkoff areas overlap. Then, we remove 7, 7', and all intervals ./ e ,c/ 
in between 7 and I'. Finally, we add the smallest interval K into ^ which contains 7 and / (i.e., 
K= (min(fl,fl').max(/7, /?'))). We proceed in this way until we obtain intervals 7j , . . . ,7^, such that 
the corresponding Obreshkoff areas do not overlap. It remains to show that the length of each of 
the so-obtained intervals is bounded by Wmin • n^+^iog"; At any stage, an interval / obtained in 
the above iteration, covers a certain number s of intervals from 7i , . . . ,7^. By induction on s, we 
prove that 

wiJ) < VVn.in • /«21og.v^ ^^^^^ ^(y) < ^^.^ . „21og«+4 ^ 

If i = 1, then trivially J =Ii for some / which proves the claim for s = 1 . An interval J covering 
5-1-1 intervals from 7i , . . . ,7^ is obtained by merging two intervals 7 and 7' which cover s\ and ^2 
intervals, respectively, where si +S2 <s+l and, w.l.o.g., 1 < < S2- From Lemma 3, it follows 
that the distance between 7 and 7' is bounded by • min(w(7'), w(7)) because the corresponding 
Obreshkoff areas overlap. Hence, J has width 

w{J) < w{I) +w{l')+4n^ ■ min(w(7'), w(7)) < w^n • {sW^°^'' + {4n^ + l)4n2iog.i ) 
If 2si > S2, then 

4„2iog.2+8^yiog2., < (84+4)n2i°g2.i < (ii+i2)Vi°s(*i+^2). 
Otherwise, we have 

and thus (3.8) follows. 

Hence, we count at most k sign variations for the intervals 7;, in total. Now, the same argument 
as used in Section 2 to show that there exists at most n' = var(/, [— 2^+ ^ , 2^+ ' ] ) special nodes also 
shows that the number of all special nodes J with w{J) < cr is bounded by k. Namely, we start 
with special nodes h,...,Im with YULi ^^^{fJi) < k, and whenever a special node is subdivided, 
the value /J := L/ var(/,7) — #{7 : var(/,7) > 2} decreases by at least one. For the result on the 
width of the intervals Ij, we consider integers ki,...,ks, with ki < k2 < ■ ■ ■ < ks, such that 

C7(zi) = • • • = (7(z„-yfcJ > 0{Zn-k,+l) = ■■■ = C7(z„-,t,_, ) > a{Zn-ki+l) = ■■■ = C7(z„). 

There exist at most kj special nodes /,■ with vi'(/,) < *''' n~4~2iog«^ ^^^^ there cannot ex- 
ist a special node 7,- with w{Ji) < Of = C7(z„). In addition, the number of special nodes with 
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w{Ji) < 0'(z„_<;| )n^^^^'°S"/4 is bounded by ki. Then, by induction, it follows that w(7,) < 
(7(z,)n^''^^^'°S"/4 because the 7,- are ordered with respect to their length. 

For the bit complexity of the endpoints of an interval / = {a,b) S T{i) n {7,}, we remark that, 
due to our construction, a and b are both dyadic numbers with modulus bounded by 2^+' . In ad- 
dition, the denominator of a and in its dyadic representation is bounded by max ( 1 , log w (/) ~ ^ ) . 
Hence, we can represent both endpoints of / with 0{t + log (7,~ ' + log^ n) many bits. □ 

We will now derive our final result on the bit complexity of Dsc^: In each step of the 
algorithm, we have to compute the polynomial fi{x) = (x+ 1)" ■ f{{ax + b)/{x + 1)), where 
I = {a,b) is the interval that is actually processed. The latter computation decomposes into com- 
puting // := f{a + {b — a)x), reversing the coefficients, and then applying a Taylor shift by 
1 (i.e., x>-^ x+l). For the computation of //, we first shift / by a, and then scale by a fac- 
tor ^— a = w(/) which is a power of two. Using asymptotically fast Taylor shift [11,32], the 
shift X 1-^- X -|- a (i.e., the computation of f{x + a)) demands for d{n^{T + logw{I)~^)) bit oper- 
ations. The scaling x>-^ {b — a) -x is achieved by just shifting the i-th coefficient of f{x + a) 
by i ■log{b — a)~^ = j-logw(/)~^ many bits. Then, the resulting polynomial // has coeffi- 
cients of modulus 2"^"^'', and the denominators of their dyadic representations are bounded by 
reversing the coefficients of //, and then applying a Taylor shift by 1, 
demands for 0(n^(T-|-logw(/)"')) bit operations. In sunnmary, the costs for computing // are 
bounded by 0(n^(T-|-logw(/)~^)). The same bound further applies to the computation of Ai and 
X2 in Step 3 (b), (3.2) of our algorithm because, in this step, we have to evaluate a polynomial of 
degree n and bitsize T at a ( T -|- log w (/) " ' ) -bit number. If 7 is non- terminal, then we also have to 
compute the number of sign variations for the intervals Bi, B2, l[ and I2. The same argument as 
above also shows that we can do so using 0(n^(T-|-logA'f-|-logw(/)~^)) bit operations. 

For an interval I e T{i) U7„ we have 

logiV/-hlogw(7)-i = 0{logNi + logw{Ij)-'^) = 0{T + loga{zi)~^ +log'^n)y, 

see Lemma 1 and 6. Thus, the total costs for all computations at I are bounded by 0(n^(T -|- 
a(z,)^')) bit operations. It remains to consider a terminal interval I which is one of the two chil- 
dren of a special node 7,. In this case, we only have to bound the cost for the computation of // be- 
cause var(/,/) < 1. Since w(/) =h'(7, )/2, the latter computation needs (5(«^(T + log w(/y )"')) = 
d{n^{'V + loga{zi)~^)) bit operations as well. 

In Theorem 4, we have shown that |7; | = 0(logn -l-logT) for all i, and thus the total costs for 
isolating the real roots of / are bounded by 

«' 

0(logn + logT) • £0(n2(T + logC7(zO"^)) = d{nh) 

i=l 

since Y.'iLi log cfe)"^ = 0{nT + Y!i=i log (7(z,)"^) = (5(nT) according to Lemma 19 in [28]. We 
fix this result: 

Theorem 7. For a polynomial / of degree n with integer coefficients of modulus less than 2^, 
DSC^ isolates the real roots of / using no more than (5(«^t) bit operations. 
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4. Conclusion 



We introduced the first subdivision method to isolate the real roots of a polynomial which 
achieves the record bound 0{n^x) for the bit complexity of this fundamental problem. In com- 
parison to the asymptotically fast algorithms by Pan and Schonhage from the 80tes which com- 
pute all complex roots, the new approach is much simpler and can be considered practical to an 
extremely high degree. The algorithm is based on a novel subdivision technique which combines 
Descartes' Rule of Signs and Newton iteration. As a consequence, our algorithm shows quadratic 
convergence towards the roots in most steps. 

So far, the approach only appUes to polynomials with integer (or rational) coefficients. In [27], 
we showed how to modify a subdivision algorithm which uses exact computation in each step 
such that it also applies to polynomials / with arbitrary real coefficients that can be approxi- 
mated to any specified error bound (bitstream coefficients). Following this approach, it seems 
reasonable to express the bit complexity in terms of the geometry of the roots zi,...,Zn of /. 
For instance, the modified version of the Descartes method from [28] isolates all real roots using 
d{n{nr + 1,)^) bit operations, where F constitutes a bound on the logarithm of the modulus of 
the roots and E : = YJi= i log c (z; ) " ' • Due to the quadratic convergence achieved by the new algo- 
rithm, we expect a corresponding bitstream version to perform the same task by a considerably 
lower number of bit operations. 
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5. Appendix 



Algorithm 1 Dsc^ 



Require: polynomial / = Y.o<i<n'^i-'(' £ with integer coefficients a,-, |a, | < 2^ for all i. 

Ensure: returns a list ^ of disjoint isolating intervals for all real roots of / 
7o:= (-2^+1,2^+1); Ar/„:=4 

.s/:={ (/o, A^/q) }; ^ := {list of active and isolating intervals} 

repeat 

(/,A^/) some element in with 7 = {a,b); delete (/,//) from jz/ 

v:=var(/./) 

if V = then 

do nothing {/ contains no root} 

else if V = 1 then 

add I to ^ {I isolates a real root} 

else if V > 1 then 

Bi:=(«,«+'f );B2:=(/,-f),Z,) 
if var(/,fii) = V or var(/,fi2) = v then 

for the unique i G {1,2} with var(/,B,) = v, add (BnA^s,) := {Bi,Nf) to ^ 

{B, contains all real roots within I } 

else 

for J = 1,2: 

fc,- :=min(max(L4Ar, • |E^J ,2),4iV, - 2); /; := {a + ik - 2) • ^,a + (A;,- + 2) • ^); 
{wj: : = a + A: • w subdivision point when decomposing I into 4Ni equally 

sized intervals; mt- is closest to Aj-; has width and is centered at nik^.} 

if var(/,7j) = v or var(/,72) = v then 

choose the smallest i G {1,2} with var(/,//) = v and add (/■,A',p := {I'i,Nf) to ^ 

{Ifvai{f,I-) = V, then // contains all roots within 1} 

else 

add (/i,A^/,) := ((a,m(/)),max(4, VA^)), (/2,Af/J := ((m(/),^),max(4, y^)) to^ 

if/(m(/)) ^Othen 

add [m (/),«(/)] to ^ 
end if 
end if 
end if 
end if 
until ^ is empty 
return 
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